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INTRODUCTION 


With  the  expansion  of  testing  conducted  under  project  BIFI,  it  is 
clear  that  a high  degree  of  automation  is  required  in  the  data  acquisi- 
tion and  reduction.  When  a Hewlett-Packard  digital  scanning  system 
became  available  for  use  in  October  1968  it  became  possible  to  automate 
to  a high  degree  the  daily  propagation  tests  conducted  over  the  BIFI 
range  (described  in  reference  a).  This  memorandum  deals  with  the  BIFI 
data  acqulsion  system  located  in  Building  36,  which  is  used  in  this 
testing  and  the  two  FORTRAN  V programs,  S1298  and  S1357,  which  are  used 
to  reduce  the  data. r-  — 


IATA  ACQUISITION  SYSTEM 


The  daily  propagation  tests  are  conducted  over  the  BIFI  range, 
shown  in  Figure  1.  The  range  has  a length  of  approximately  19  nautical 
miles  and  has  a depth  of  about  120  feet  through  most  of  its  extent.  At 
Block  Island,  a projector  tuned  at  1702  Hz  is  bottom  mounted  at  a 
55  foot  depth  at.  point  S.  The  receiving  hydrophone  is  bottom  mounted 
near  Fishers  Island  in  155  foot  of  water  at  point  H.  The  receiving 
station  at  Fishers  Island  is  connected  by  means  of  data  transmission 
lines  to  the  Data  Acquisition  and  Reduction  Center  in  Building  36  at  the 
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laboratory,  where  the  signals  are  recorded  and  analyzed.  A reference 
signal  corresponding  to  the  signal  transmitted  at  Block  Island  is 
received  at  the  Laboratory  from  the  Block  Island  station  by  means  of 
the  Block  Island  telephone  data  line.  Several  frequency  sensitive 
reed  relays  are  connected  in  the  receiving  circuits  at  Fishers  Island 
and  these  permit  remote  control  and  calibration  of  the  receiving  system 
from  the  Laboratory  via  the  data  transmission  lines. 

The  receiving  system  in  Building  36  is  shown  in  Figure  2.  The 
signal  from  the  receiving  hydrophone  is  received  via  the  Fishers  Island 
telephone  data  line,  the  CW  signal  is  passed  through  either  a 2 Hz  or 
0.1  Hz  bandwidth  filter  prior  to  digitizing.  The  2 Hz  bandwidth  filter 
is  used  in  the  winter,  spring,  and  fall  when  the  signal  to  noise  ratio 
of  the  received  signal  is  large  and  there  is  a large  amount  of  fre- 
quency smear  in  the  signal.  The  0.1  Hz  bandwidth  filter  will  be  used 
in  the  summer  months  when  the  signal  to  noise  ratio  of  the  received 
signal  is  smaller  than  in  other  seasons  and  due  to  a lack  of  dependence 
of  the  received  signal  upon  the  sea  surface,  thus,  there  is  very  little 
frequency  smear.  The  envelope  of  the  filtered  signal  is  then  sampled 
by  the  Hewlett-Packard  scanner  and  digital  voltmeter  system  and  the 
digital  value  punched  in  IBM  cards  at  a rate  of  approximately  3.5 
samples  per  second.  The  entire  receiving  system  has  a dynamic  range  of 
better  than  50  db.  Over  a range  greater  than  40  db  the  system  is  linear 
to  at  least  .2  db  per  10  db  step.  The  system  noise  is  low  enough  to 
permit  measurement  of  signals  of  as  low  as  -60  db  relative  to  a 1 Dyne 
per  square  em. 

During  the  daily  propagation  tests,  pulses  45  seconds  long  are 
transmitted  from  Block  Island.  During  each  one  minute  interval  the 
signal  is  on  for  about  45  seconds  and  off  for  approximately  15  seconds. 
The  recording  sequence  during  the  tests  is  shown  in  Figure  3.  First 
the  Block  Island  reference  signal  is  recorded  on  IBM  cards,  prior  to 
the  first  signal  sequence  being  received  at  Fishers  Island,  for  more 
than  one  minute  or  one  cycle.  Then  the  received  signal  at  Fishers 
Island  is  recorded  for  15  minutes  on  cards.  Finally,  a calibration  of 
the  system  is  performed.  The  calibration  level  is  changed  in  equal 
increments  and  a seme  number  of  values  are  record jd  at  each  level.  The 
format  of  the  data  on  the  cards  is  four  digits  representing  an  integer 
followed  by  one  digit  which  represents  a negative  exponent  of  the  number 
10  by  which  the  first  number  is  to  be  multiplied.  Next  the  data  is 
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PROGRAM  S1298 

USL  Program  SI 298 , written  in  FORTRAN  V language,  is  designed  to 
determine  propagation  loss  and  variance  of  signals  received  during  the 
daily  propagation  runs.  The  data  obtained  is  arranged  as  follows: 

1.  One  input  data  card. 

2.  Calibration  cards  with  higher  calibration  levels  preceding 
lower  levels. 

3.  The  signal  data  cards  in  the  order  in  which  the  data  is  taken, 
as  described  previously. 

The  format  of  the  input  data  card  and  the  quantities  inputted  are 
shown  in  Table  1.  As  seen  in  Figure  3 the  first  data  reoorded  consists 
of  the  Block  Island  reference  signal.  This  signal  is  used  in  program 
S1298  for  two  purposes.  First,  since  a complete  signal  cycle  is  one 
minute  in  length  from  point  A to  B,  we  can  compute  the  sampling  rate 
per  minute  as  the  number  of  samples  in  the  period,  EMIN.  The  reference 
signal  is  also  used  to  establish  the  time  at  which  the  signal  is  first 
transmitted  from  Block  Island  (point  B).  From  the  known  propagation 
time  (input  parameter  PROT  which  is  24  seconds  for  the  present  BIFI 
range)  and  the  time  after  the  beginning  of  the  received  pulse  at  which 
analysis  is  to  commence  (input  parameter  ENC  which  is  5 seconds  for 
present  analysis)  the  program  determines  point  C and  the  signal  is 
analyzed  from  point  C to  point  D (DUR  seconds  long  which  is  30  seconds 
for  the  daily  propagation  runs).  The  ambient  noise  is  analyzed 
commencing  at  point  E.  The  distance  from  point  B to  point  E is  the 
input  parameter  RAMB  (74  seconds  in  the  present  analysis).  The  analysis 
of  ambient  is  conducted  from  point  E to  F (DURA  seconds  long  which  is 
5 seconds  in  the  present  case).  Analysis  of  succeeding  pulses  is 
accomplished  by  incrementing  points  C and  E by  one  minute  (EMIN).  The 
, number  of  pulses  to  be  analyzed  is  given  by  INUM. 

For  each  pulse  analyzed, calculations  are  performed  to  measure  the 
maximum,  minimum  and  average  values  of  voltages  over  the  pulse,  and  the 
variance  of  the  log  of  the  voltage  over  the  pulse.  Also  calculated  is 
the  average  value  of  the  voltage  representing  the  background  noise. 

A number  of  calibration  levels  (input  parameter  NCAL  is  the  number 
of  levels)  are  analyzed.  Each  level  has  CAL  measured  values  and  the  total 
number  of  calibration  values  analyzed  equals  ICAL.  The  calibration 
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equivalent  level  in  db  relative  to  1 DYNE  per  square  cm  of  the  highest 
calibration  level  is  given  by  BQIJL  and  each  level  is  decreased  by  ECDP 
db  in  each  of  the  other  calibration  levels.  The  calibration  equivalent 
loss  of  the  highest  calibration  level  is  given  by  EQUV, 

EQUV  = N,  - EQUL 


where  Na  is  the  source  level  and  EQUL  is  the  calibration  equivalent 
level,  both  in  db  relative  to  1 DYNE  per  square  cm.  Each  loss  level 
is  incremented  by  ECDP  db  in  each  of  the  other  calibration  levels. 

For  each  calibration  level  an  average  voltage  level  is  calculated  and  an 
array  SIGJ(JJ)  is  set  up  such  that 

SIGJ(JJ)  = 20  Log  ( V ) 

where  V is  an  average  voltage  and  JJ  = (1...NCAL)  where  SIGJ(l)  is 
the  highest  calibration  and  each  succeeding  value  at  the  array  repre- 
sents a level  with  a calibration  equivalent  level  of 

EQUL  - (J./-1)  * ECDP 
or  a calibration  equivalent  loss  of 

EQUV  + (JJ-1)  * ECDP 

Having  calculated  the  calibration  levels  and  the  maximum,  minimum, 
and  average  of  each  pulse  received  level  it  is  possible  to  calculate 
the  propagation  loss  associated  with  the  later  three  quantities.  For 
a given  level  of  one  of  these  quantities  propagation  loss  is  calculated 
by  either  interpolating  between  or  extrapolating  beyond  calibration 
levels  associated  with  distinct  calibration  levels  equivalent  losses. 

Let  x be  the  quantity  whose  associated  propagation  loss  Px  is  desired  to 
be  measured. 

If  x lies  between  two  calibration  levels  then  we  interpolate 
between  those  values  and 

Px  = EQUL  + ECDP  • MM  - (10.0  • LOG  (X)  - SIGJ  (LL)  • EDCP/ 

(SIGJ(MM)  - SIGJ(LL) ) 

where  SIGJ  (MM)  and  SIGJ  (LL)  are  the  values  of  SIGJ  corresponding  to 
calibration  values  above  and  below  x respectively. 
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If  x Is  greater  than  any  calibration  level  then  we  extrapolate 
beyond  the  upper  calibration  value  and 

Px  = EQOL  - (10.0  • LOG(X)  - SIGJ(l))  • ECDP/  (SIGJ(l)  - SIGJ(2)) 

If  x is  less  than  any  calibration  level  then  we  extrapolate  below 
the  lowest  calibration  value  and 

Px  = EQOL  + (NCAL-1)*ECDP-(10.0*L0G(x)-SIGJ(NCAL) )*ECDP/SIGJ 

(NCAL-l)-S'GJ(NCAL) 

If  the  number  of  calibration  values  NCAL,  is  one  it  is  assumed  that 
the  recording  system  is  linear  between  x and  the  calibration  level  and 

Px  = EQUL-10.0*L0G(x)+5IGJ(l) 

If  x is  the  average  background  noise  level,  P_  is  the  propagation 
loss  which  would  be  measured  for  a signal  with  a level  equal  to  the 
background  noise.  This  quantity  AMB2 , is  a measure  of  the  mariwim 
propagation  loss  which  can  be  measured  for  a given  level  of  background 
noise.  The  background  noise  level  AMB  is  related  to  AMB2  by  the 
relationship 

AMB  = - AMB2  + EQUL  + BQOV 

If  the  maximum  or  minimum  level  of  signal  plus  noise  falls  below 
the  background  noise  the  quantity  is  set  equal  to  200  db  propagation 
loss.  If  the  average  value  of  signal  plus  noise  does  not  exceed  the 
background  noise  then  the  average  value  of  propagation  loss  is  set 
equal  to  AMB2  and  GTR  is  set  equal  to  one.  This  effectively  states  that 
the  average  propagation  loss  is  greater  than  AMB2. 

For  each  pulse  the  maximum,  minimum,  and  average  values  of  propaga- 
tion loss  in  db  are  printed  out.  Also  printed  out  is  the  variance  of  the 
propagation  loss  in  db  over  the  pulse  length  as  well  as  the  ambient  in 
db  relative  to  1 dyne  per  square  cm,  and  AMB2  is  db. 

Also  printed  out  are  the  average,  of  the  number  of  pulses  analyzed 
of  the  maximum,  and  average  propagation  losses  and  the  average  variance 
over  a pulse.  The  variance  of  the  maximum  and  average  propagation  loss 
measurements  of  the  number  of  pulses  analyzed  is  also  printed  out.  These 
quantities  are  also  punched  on  one  output  card,  the  format  of  which  is 
shown  in  Table  2.  Also  punched  out  on  this  card  are  sea  state,  and 
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wind  speed  at  the  time  of  the  tests  and  the  value  of  GTR.  This  card 
can  be  used  as  data  input  to  Progran  S1357. 

There  are  two  printed  outputs  in  program  S1298  which  can  be  used 
as  a check  of  the  accuracy  of  the  measurements  taken.  First  the 
voltage  level  in  db  of  each  calibration  level  is  printed  out  ao  that 
the  linerity  of  the  data  acquisition  and  reduction  system  can  be 
verified.  The  second  check  concerns  the  exponents  of  all  input  data. 
Due  to  the  nature  of  the  Hewlett  Packard  system  the  exponents  of  all 
data,  taken  with  a given  maximum  input  level  setting,  have  the  same 
value  for  the  exponent.  This  value  is  inputted  as  the  quantity  EXP. 

If  the  value  of  the  exponent  of  any  data  point  differs  from  E3P,  then 
the  value  of  the  exponent  and  its  location  in  the  data  deck  is  printed 
out.  Thus  it  is  insured  that  a faulty  alignment  of  data  card  columns 
in  the  data  can  be  detected. 


PROGRAM  S1357 

USL  Program  S1357,  written  in  FORTRAN  V language,  is  designed  to 
measure  averages  of  the  propagation  loss  and  variance  data  obtained  from 
program  S1296  and  to  plot  these  as  well  as  theoretical  predictions  when 
applicable  on  the  cal comp  plotter  in  a form  such  that  the  graphs  may 
be  used  in  technical  reports.  The  data  is  arranged  as  follows: 

1.  One  input  data  card. 

2.  Input  cards  containing  data  of  theoretical  predictions. 

3.  Output  cards  from  Program  S1298. 

4.  Cards  containing  a list  of  sea  states  to  be  used  as  input 
parameters. 

5.  Cards  containing  a list  of  wind  speeds  to  be  used  as  input 
parameters. 

The  format  of  these  cards  is  shown  in  Tables  II  to  V. 

Program  S1357  is  designed  to  generate  as  many  as  ten  Calcomp  plots 
from  a set  of  data  generated  by  many  propagation  tests  whose  data  has 
been  analyzed  by  means  of  program  S1298. 
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A sample  plot  of  one  of  the  ten  graphs  which  can  be  obtained 
Is  shown  In  Fig.  4.  All  values  of  propagation  loss  are  plotted 
versus  sea  state  and  wind  speed.  Average  values  of  propagation 
loss  are  also  plotted  against  sea  state  and  wind  speed.  The  values 
are  compared  with  those  predicted  by  the  COLOSSUS  equations.  In  the 
second  set  of  curves  standerd  deviations  of  each  average  value  are 
plotted.  A special  symbol  is  generated  for  data  whose  value  for 
GRT  is  one.  This  symbol  indicates  that  the  value  plotted  is  an 
indication  of  a minimum  propagation  loss  which  could  be  measured. 

Average  propagation  loss  determined  by  measuring  an  average 
pressure  level  in  a single  pulse  is  compared  to  that  determined  by 
measuring  a maximum  pressure  level  in  a pulse  in  two  graphs.  Both 
are  plotted  versus  sea  state  and  wind  speed. 

The  average  of  the  variances  of  propagation  loss  measured  over 
individual  pulses  is  plotted  versus  sea  state  and  wind  speed.  Values 
for  which  GTR  equals  one  are  not  plotted. 

The  variance  between  pulses  of  measured  propagation  loss  in  a 
given  run  determined  by  measuring  an  average  pressure  level  in  a 
single  pulse  is  compared  with  that  of  propagation  loss  determined 
by  measuring  maximum  pressure  level  in  a pulse.  This  is  done  in  two 
graphs  where  both  are  plotted  versus  sea  state  and  wind  speed. 

Values  for  which  GTR  equals  one  are  not  plotted. 

If  ISZ  (see  Table  III)  is  set  equal  to  one,  only  the  four  graphs 
in  which  propagation  loss  is  determined  from  average  pressure  levels 
are  plotted.  • If  ISZ  equals  zero,  all  ten  graphs  are  plotted.  If 
ISEA  is  set  equal  to  one,  labels  pertinent  to  "winter  data"  (small 
positive  to  weak  negative  velocity  gradients)  are  put  on  the  graphs 
while  an  ISEA  of  zero  corresponds  to  "summer  data"  (weak  negative 
to  large  negative  velocity  gradients). 

Tables  are  also  printed  out  which  list  sea  state  or  wind  speed, 
average  propagation  loss,  number  of  values  at  the  particular  wind 
speed  or  sea  state  considered  and  standard  deviation  of  these  values. 
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SHttARY 

The  current  BIFI  system  of  data  acquisition  and  reduction  of  the 
daily  propagation  tests  brings  a high  degree  of  automation  to  the 
processing  of  the  acoustic  data.  In  the  future  it  is  planned  to 
transmit  signals  at  least  in  three  different  frequencies  and  to  receive 
these  signals  with  hydrophone  arrays  located  at  Fishers  Island  and 
Watch  Hill  and  connected  via  telephone  data  links  to  the  data  reduction 
laboratory  at  DSL.  It  is  planned  to  purchase  Kennedy  1406  digital  tape 
recorder  to  replace  the  IBM  card  punch.  This  will  increase  the  sampling 
rate  attainable  from  about  15  samples  per  second  to  100  samples  per 
second.  Since  the  Hewlett  Packard  System  can  scan  through  a maximum 
of  100  channels  with  small  modifications  to  the  receiving  systems  and 
programs  one  will  be  able  to  process  data,  in  a maimer  similar  to  that 
described  in  this  report,  at  many  frequencies  and  many  receiving 
hydrophones . 


WILLIAM  G.  KANABIS 


REFERENCES 
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concerning;  DSL  Technical  Memorandum  No.  2211-18-68 
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FORMAT  (2IS,  FS.l,  6F5.0,  F5.1,  2IS,  2FS.0,  FS.l,  FJ.O) 

INPUT 


PARAMETER 

COLUMNS 

INUM 

1-5 

Number  of  pulses  to  be  analysed 

ICAL 

6-10 

Total  number  of  Cal  values  to  be  read  in 

EQUV 

11-15 

Cal  equivalent  loss  of  first  (highest  value) 
calibration  (source  level  - cal  equivalent 
level  (DB)) 

PROJ 

16-20 

Propagation  time  in  sec  between  source  and 
receiver 

ENC 

21-25 

Time  after  leading  edge  of  received  pulse  at 
which  summation  should  begin  (SEC) 

DDR 

26-30 

Length  of  signal  sample  to  be  analysed  (sec) 

CAL 

31-35 

Number  of  values  in  individual  calibrations 

RAMB 

36-40 

Time  after  leading  edge  of  transmitted  pulse  at 
which  noise  should  be  sampled  (sec) 

DURA 

41-45 

Length  of  ambient  desired  to  analyse 

EQUL 

46-50 

Cal  equivalent  level  of  first  (highest)  calibra- 
tion level  (db) 

ISIG 

51-55 

Total  number  of  values  of  signal 

NCAL 

56-60 

Number  of  calibration  levels 

ECDP 

61-65 

Increment  in  levels  of  calibrations  (db) 

SS 

66-70 

Sea  state 

WSPD 

71-75 

Wind  speed  (mph) 

EXP 

76-80 

Value  of  exponent  in  all  data 

TABLE  I 

INPUT  CARD  FOR  PROGRAM  S1298 
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FORMAT  80  (rio.5) 


PARAMETER 

COLUMNS 

- 

SS 

1-10 

Sea  State 

WSPD 

11-20 

Wind  Speed 

AVE 

21-30 

Average  propagation  lose  measured  over  pulses 

EMAX2 

31-40 

Average  of  maximum  propagation  loss  measured 

AVAR 

41-50 

Average  variances  over  pulses 

VMAX 

51-60 

Variance  of  maximum  propagation  loss  between 
pulses 

VAV 

61-70 

Variance  of  average  propagation  loss  between 
pulses 

GTR 

71-80 

l=propagation  loss  greater  than  AVE 

O=propagatlon  loss  equal  to  AVE 

TABLE  II 

OUTPUT  CARD  FROM  PROGRAM  S1298 


USL  Tech  Memo 
No.  2211-91-69 


FORMAT  315,  F5.0,  315 
INPUT 


PARAMETER 

COLUMN 

NC 

1-5 

Number  of  prediction  values 

ND 

6-10 

Number  of  data  cards  from  program  S1298 

NS 

11-15 

Number  of  sea  states  inclusive  between  lower 
and  highest  values  in  increments  0,  0.5,  1,2, 
3 

EST 

16-20 

Propagation  loss  value  at  origin  of  graphs 

NV 

21-25 

Number  of  wind  speeds  Inclusive  between 
lowest  and  highest  values  in  increments 

0,  2.5,  5,  10,  15  

ISZ 

26-30 

Number  of  graphs  0 = 10  1=4 

ISEA 

31-35 

Season  data  was  taken  0 = summer;  1 = winter 

TABLE  III 

INPUT  CARD  PROGRAM  S1357 
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FORMAT  5F10.5 

NC  = NUMBER  OF  CARDS 


INPUT 

PARAMETER 

CSS 

CWS 

COLU 

COLL 

CW2 


COLUMNS 

I- 10  Sea  State 

II- 20  Wind  Speed 

21-30  Upper  Propagation  loss  prediction 

31-40  Lower  Propagation  loss  prediction 

41-50  Wind  Speed  (2) 


TABLE  IV 

INPUT  CARDS  WITH  THEORETICAL  PREDICTIONS 
PROGRAM  S1357 
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NS  = NUMBER  OF  SEA  STATE  VALUES 
NW  = NUMBER  OF  VIND  SPEED  VALUES 
FORMAT  8F10.5 
INPUT 

PARAMETERS  COLUMNS 

Z or  ZZ  1-10  ZZ  is  wind  speed  (mph) 

11-20  They  are  inputed  in  ascending  order  of 
magnitude 

21-30 

31-40 

41-50 

51-60 

61-70 

71-80 


TABLE  V 

INPUT  CARD  LISTING  SEA  STATES  OR  WIND  SPEEDS 
PROGRAM  S1357 
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luttl  ulTY  CALCUi-mI  I Oil  FoR  LONo  PUL-»ES  S\arW 

Ui.  ri  bK'N  Lio  ( 10o0ij>  . SIGC(IUOO)  • ESIGl  10000)  ,ESlGc  ( 1000 ) 

1 , F ( 10000 ) , Slr7J<l()0) 

WRIT*  (4,19) 

ih  i Ji\i<IA  f ( lUtiH  MAX#  till##  AvEk  are  PhOP  LObS  IN  Dd.  AMH  IN  OB  REL  DYNE 
1 PlR  JsoUAkE  CM.  VAR  Ii,  VARIANCE  OF  LOG  PRES  SQUARED) 

RL>.D  ( J,  1)  INUN  » II AL, lQuV,PROT,F.Nc»UUR# CAL, RArtd, DURA, EQUL.ISlG 
1 » KoAL  . Etop  » S$  » MSPD  » EXP 

wRaTE (4#1)  IhUM,  ICAL#LQUV,PROT»f.NC»UUR#CAL»RAMB,OURA#EQUL,  ISlG 
1 # NLmL  # ECi  P » SS  » WSPD  » EXP 

1 Format  ( <.  1 b > Fb • 1 # Op S. u » Fb * 1 # 2 1 b# 2F‘j .0,F5.1,F5.0) 

Rt».D  (o»2)  ( bluC(J)  , lSIOC(J),  j=1#ICAL  ) 

2 F Oh VAT  Clb(F4.G,Fl,0)  ) 

KumU  (3,2)  ( SI  v»  ( J ) » Lb  I G ( J ) # j=  1.  ISlG  ) 

Do  70  J = 1#  ICAL 

IF  (EblGC(J)  - 3.0  ) 72,70,72 

72  WRITE  (4# 7b)  J . LSIoC(J) 

75  Format  ( lb#  f5.c  ) 

70  CONTINUE 

DO  71  J = 1#  ISlG 

If  ( ESIG(J)  - EXP  ) 73»7l#7 

73  WRITE  ( 4 » 7 b ) J » ESlG(j)  „ n,..T  rTY  PRACTICAL 

7b  FORMAT  <bX»  15#  Fb, 0 ) THIS  PAGE  IS  BEST  ^ALIT*  , 

71  CONTINUE  rawiCOPYEUhlUSHiSOTO^0  

Uo  60  Jj  = i » NCAL 

60  SioJ(JJ)  = 0,0 

ICAL  = ICAl/uCAL 

DO  61  JJ  =1#  NCAl 

DO  3 J = 1#  ICaL 
JL  = J + ILAl*CJo-1) 

3 SluJ(JJ)=bIoJ(Ju)-KSIt.Cr  JL)/10.*n£bIGC(JL)  )**2.0 
SlOJ(JJ)  = 10.0*AL0olu(bIGJ(JJ)/CAL) 

ib  Format (Flu.b) 

bl  Wr1jE(4, lb)bloJ(Jj) 

DO  4 J=  1 # 30  0 

F (o+l ) = SIG(J+1)  / I0.0+*  EbIG ( J+l ) 

30  format  (Fio#b) 

Fiw)  = blG(J)/i0,0+*  EblG(J) 

IF  ( FCJ+i)  -6 ( o ) -.2  ) 4,4,5 

4 Coi.TINUL 

5 IblR=  O + 1 

Do  o K — 1 , 3 0 u 

J = ibTR  + h 

F(w)  = blolJ)/10.P**  FblG(J) 

F (u  + 1 ) = Si G ( o + l ) / lr..U**  ESIG  ( J+l ) 

If  ( F(  J)  - F ( J +1)-  .2  ) 6,6,7 

6 Co.  TINUL 

7 IblP  - + o 

DO  H r*= 1,  30u 

o = iS  TP  ♦ r 

FVo)  = r»I6(J)/l0.0  + * L,lo(J) 

F ( „ +i ) = blb(w+l)  / 10.0**  LSIG ( J+ l ) 

If  ( M w +1)-  r ( Jl-,2)  6.6,9 

r C 01 . T I UoE 
9 Ijlti  = +0+1 

WRiTt  (4,1b)  laTh  »IbTp  ,IbTH 

i<>  Forma  f (3110)  USL  Tech  Memo 

1m  * J = lb  1 6 * I-»Tt<  No.  2211-91-69 

A— 1 


6»6»7 


SI6(J)/I0.0+*  L,lo(J) 
■jlbCw+l)  / lO,  j*+  LSIGiJ+l) 


J)-#2)  6,6,9 
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I I fij,  vl  ( boi il 


_ Li-.ii**  (I'i'O  I ♦ ii<(.  ) /o>JiO 
z i(mi“U*i  !•  11  ,/hu  • u 
- bUr\  • ./uli  . U 

> z JUhA*i.MIi;/bD.lj 

L-U.i  t If.So  »I00R 

i (Silo) 


f-*  1 1 1 


, liJUHA 

AVER 


VAf-’ 


AK« 


AW1  Z 

11 . 0 

Lei.  X. 

Z u.U 

«(i  1 

z 0..1 

» ’) 

z U . 0 

V.„.A 

C 

• 

c 

V/«\> 

z 11 . 0 

Uu  1(. 

Kz  1 1 lulu 

Lfi*>  * 

=u  .0 

tf-iiM 

Zj  UUOlUlU  . V 

To  1 b 

-u  . 0 
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T o i — 1» « 0 

T Uc.  - 0 

TjIa  =o.O 

Ou  1 1 L—  1 » *t-  Oh 

>1  = L ■*  (iv-1)  *1-  'IU 

J z 1 jTl>  * .SSo  + K' 

Tkw.)  z S10(J)/10,0»*  E'SIO(J) 

lh  ( ) (J)**k!.u  ,01  • U-V'X ) LNiAX  = ( F(J  ))**2.0 

Ilr  1 MJ)**£.o  .LT.LMN)EMlM  = ( F(J  ))**2.0 

rwl  = 10T  ♦ ( M J ) ) **2.0 

To..  z ToP  ♦ cO.O*Al.O(,1o(K(  J)  ) 

I O I b=  roTj+  (co.  0*ALoGl0  (F  ( J)  ))*♦«;.  o 

11  Cut  rl'»bc. 

Iaj  IP  |4-1*  AiX'RA 

liZ  M ♦ (K-l  ) ♦I.  ili- 

J - lb  To  ♦ ll-AMb  4 t) 

f [s.  ) - SlMjl/lO.O**  E.SIo(J) 

To  I A=  7 o T A + ( ► (J  ))**2,0 

12  COxTlhUt 
Uo.  zlu'm 
Uor  .(Z  lolUivl 

tooL  = UCAL  - ) 

lr  ( uCAl  **  1 ) 

07  A.ij  z r.oOo  - 1 0 . 0*Al_0ol  0 ( TOTA/l’URA ) SlGj(l) 

Go  T O io 
S»/»  l.<o  lL  - ltNlM 

lr  ( 1U  ,0*AlU01()  ( IUTA/i^OkA  ) , uT  .b  IG  J (LL  ) ) o0  To  4b 
o s Co*  tl* is»- 

h iu  2 z fc. Goo  *uC/  L^fc.C  J»  - ( 1 0 . U*m1.00  1 1)  ( ToT  A/iMiHa 

l+_k  i»/('.,IuJ|U/R-l)  - SlGj(  1C«|.)  ) 

Go  I M. 

'♦S  • >0  'U  Z i«L.»U  » 1 r — 1 

lr  1 lU.il*.«o'J01i-(l  J|A/.  I|  A j ,l.T,'  lGo(f  “)  ) :>J  I J iit, 

'♦/’  f I 

- * ^uu  -il  # u ♦ /»t  j,  u ( I jl/i/!  •’  ) -Sicj( ) ) ) ST  » J(  1 ) 

i " i i L w ( * ■ ) ) 

ri 

*•  /’■ . 1 » h.  ^ L w t h • 1 4 1 ( T 01  A ) ••  3 J >»*  ( I L 

»)>•  n.-iV  ( ol - o ( • ' > - RiidllU)  USL  Tech  Memo 

A-2  Mo.  2211-91-69 
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ct*  IF  (LI./.X  - ToWUUHA  ) 20,20,21 

cO  LmmX  S +cOU»U 
go  To  22 

l1  IF  ( uCAL  - 1 > 91,91*92 

Vl  tN„X=  LOUL  - 10  « U* ALOo 1 u ( LMAX-  TOTA/DURA  ) + SIoJ(l) 

Go  To  22 

92  Go  32  LL  - 1 * MLaL 

If  (lU.O*ALOblO(tN.AX)  ,GT.SIGJ(LI  ) ) GO  TO  33 

32  Cu.-rlNUL 

Em„x  = FOUL  +uCML*tC0H-( 10.0*ALOG1U(EMAX-ToTA/DURA)-SIGJ<NCaL) )* 
lEuuP/ (bIGj(UCAL  -1)  - SlGj(NCAL)  ) 

Oo  TO  22 

33  L»0  34  mm  = NCAL  ,1,-1 

IF  ( lU.O*ALOG!0(tMAX) .LT.blGJ(KM) ) GO  TO  35 

34  CONTINUE 

EMmX  = EOUL  -(lo,u*ALoGlO(LMAX  -TOTA/UURA)  -SlGJ ( 1 ) ) *ECDP/ ( SlGJ ( 1 ) 
1 - S1GJ(2) ) 

Go  TO  22 

3b  ErtAX  = FOUL  + ECDF'  + mM  - ( 10.0*ALOG10  (EMAX-TOTA/DURA)  -SlGJtLL 
1))*  ECUP/  (blGj(MM)  - SIGJ(LL)) 

22  IF  (EMIN  - TOTa/OUHA  ) 23*23,24 

23  EMN  = +200.0 

Go  TO  2b 

^4  IF  ( NCAl  - 1 ) 93,93*94 

93  Emins  FOUL  - 10.0*ALOGl0<EMIN-  TOTA/DURA  ) F SIGJ(1) 

GO  TO  2b 

94  UO  37  LL  = 1 , NCAL 

If  (10.0*ALOG10(EM1N) .GT.SIGJ(LL) ) GO  TO  43 
37  Cui.TlNUt 

Emin  = EOUL  +FCAL*LCDp-< 10.0*ALOG10(EMIN-TOTA/UURA)-SIGJ<NCAL) )* 
1Ecup/(S1GJ(NCAL  -i)  - SjGj (NCAL ) ) 

GO  TO  cb 

43  Uo  40  mM  s NCAL  ,1,-1 

If  ( 10.0*ALOG10 (EMIN) .LT.blGJ(MM) ) GO  TO  46 
40  Continue 

EmIN  e EOUL  -Uu.O  + ALOGiO  (LMIN  -TOTA/UURA)  -SIGJ(I)  )»LCDP/(SlGJ(l) 
1 - S 1GJ ( 2 ) ) 

GO  TO  2b 

46  EMIN  = LOUL  + ELUF'*mM  - ( 1 0 . U + AL0G10  ( LMIN-TOTA/DURA > -SlGJILL 
1))+  tCUP/  (bloJ(MM)  _ SIGJ(LL) ) 
tb  IF  (IOT/OUR  - ToTA/OuRA  ) 26,26,27 

26  A v lRS  aMu2 

G ( k s 1.0 


GO  10  26 

£.1  If  ( imL Al  - 1 ) 9b,  9b  *96 

9b  AwuRS  EOUL  - 10.0*AL001ij(TOT/UUR  - lOTA/OUFA)  +5lGJ(l) 
GO  10  2u 


96  L»u  36  lL  = 1 * nC/'L 

1 r (10.li*ALOGlG(Tol/riUR)  .or.SIv^J(uL)  ) GO  TG  44 
36  Coi.T  INUL 

AVLR  s t'uUL  +LC»'L*t  CUP-  ( 1 U , 0*  aLOG  L IJ  ( TOT  / OUf'-T  OTa/OUR  A ) -SlGJ  ( Nc  AL  ) ) 
1 *LUJP/ ( bloJ ( i ICAl-1 ) - SIGj(nCaL)  ) 

Go  TO  20 

44  (Go  41  =.  NCAL  *1,-1 

It-  ( 10 ,0*AlOG10  IToT/UoR ) ,lT.SIOj  (f  ,M)  ) 00  To  4/ 

4 1 Co>  TluuL 


A vi  4 = F OOL 
1 ( b i GO ( l J — 


“ ( It , u*ALooiO  ( TOT/t.l  lfi-I  oTa/OURA  ) -bloJ(  1 ) >+ECDp/ 


SlGJ ( 2 ) ) 
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47  A^uH  - LJUL  ♦ ELUP*MM-( 10.()*Alug10(T0T/DUK-T0TA/DU«A)  -SlGJ(LL 
Dm  LLbP/  (blOJiMM)  - SIGJ(LU) 
til 1 VHl  = TOTb/uUK  - (T0d/L)UH)**2.0 
Ai^u  ♦EdUL  fLUUy 

Avt.  = mVE  ♦ AVER 
Arf^ri  - AVAM  ♦ Ai*)H2 


A„/*R  = AVAR  + VAR 

EW«X«  = EMAX2  + EM A A 
VM*X  = VMAX  + EMmX**2.0 

V A v = VA V + AvER**2.0 

WKiTL  <4»13)tMAA,EMfi  iaVER»  VAR,  AMB,AMGi> 

13  Format  (6(Fio.b) ) 


in 


Cui.TlNuL 
A vu  - 

A^aK  = 
AVAM  = 
Em«X2  = 
VMmX  = 
V*v 


ave  / H-oat(inuM) 

AVAR/  FLdAT(lNuM) 

AVAM/  FLOAT (1NUM) 

EMAX2  / FLOAT (INUM) 
VMAX/t-L0AT(INuM)  “ EMaX2**2.0 
VAV/FLOmTIINUM)  - A vfc **2  • 0 


ral8MralsB^mmcnoBi 

rROM  COPX  10  »«• 


14 


31 


fa5 


WRi TE (4  * 14 ) EMAA2  » aVl  » Avar  >avam 

FoKMAT(F1U.5»10a#2F10.5  »10X  » FlU.5  ) 

WRITE  (4»31)  VI'  AX  t \t*\j 

Format  (Fio.b, iox»fio.5> 

write  i 9»ob)  bs#  wspu  • ave  . emax2  » 

Format  ( oFio.b  ) 

£nl 


AVAR  # VMAX  » VAV  ,oTR 
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tun  plots  L>»ac~\ 

DIMENSION  LOLO(IOO)  , LOLL  ( 1 0 0 ) . Ai2000)  . Y<2000)  , oAThIIO^J  , 
1C55(100)  .CWSdUO)  , FN(tUO)  i \<(lUu)  » 21100 WSUuO ) »Sl IlOO ) # 

1 S2(l00),  W(iou)  , CM2 (100>  » XX (2000)  t 22(100)  » FM(IOU)  * 

1 YY  (2000)  . YYy  (2y00)  * YYYK2000)  , YYY2(2u00)»  W*<100).  MWM(100 


GTRUOOU) 

» XXU(3) 


»15HM 

. nw 


9o 


13 


o.o 


-3  ) 


1)  » MWWl(lOU)  . ***2(100)  , 

1 .XL(3)  . YL(3)  , 2L(3) 

cal.l  plots  (OA t a ( i > ,1024*0) 
call  plot  (O.o,  0.0  * -3  ) 

CALL  SYMbOL (5.0  , i.o  » 0.2l 
READ  (3.1)  NC  * ND  , NS  »EST 

1 FGaMaT  (315  »F5.u  .315  ) 

KtAD  (3.2)(CSSH)  . CMS(I)  .COLU(l) 

2 FOKMAT  (5Fl0.b) 

KLAD  (3,3) (X( Z ) .Xx(l).  Y(l)  . YY ( I ) 
1(1). 1=  l.ND) 

3 FOKMAT  (8F10.5) 

WRITE  (4.96)  XX(7) 

WRITE  (4.96)  Ax (8) 

FORMAT  (F10.S) 

DO  30  1=  1.2 

call  plot  (20.0. 

YCND  +1)  = EST 

Y (ND  +2)  =10.0 

C0lU(NC+1)=  EST 
COLU(NC+2)=10.0 
C0LL(NC+1)=  ESI 
COLL (NC+2) =10.0 
IF l I t Eq ,2  ) GO 
X(NO  +1) 

X (ND  +2) 

CSS(NC+1 ) 

CSS(NC+2) 

GO  TO  21 
Xa (ND  +1)  = 0.0 
XX ( ND  +2)  =5.0 
CWS(NC  +1)  = 

Cm5(nC  +2)  = 

CW2(NC  +1)  = 

Cm2 (NC  +2)  = 

CALL  LINE  (XX. ». 


. GTA(IOO) 
22L(3) 


aANABIS  S1357.90.0 
.152  . 1SEA 


. WL(3) 


» 15  ) 


. COLL(I) ,CW2(I> . 1=1. NC) 

, YYY  ( 1 ) » YYY1U).YYY2(I),GTK 


=0.0 
=1.0 
= 0.0 
= 1.0 


TO  13 


S?#** 


o.u 

5.0 
0.0 
5.0 
Nu, 

DO  71  L = l.ND 
IF  (GTK(l) .LT.0.5) 
XXL ( 1 ) = XX(L) 


l.-l.  5 ) 

GO  TO  71 


XXL ( 2 ) 
XXL (3) 
YL(1)= 


= 0.0 

= 5.0 

Y(L)  +1.0 


YL(2)  = 

LST 

YL(3)  = 

10.0 

call  line 

(XXL  . YL 

.1 

.1 

» 1 » 2j, ) 

n 

CONTINUE 

CALL  LINE 

(CWS  , CoLL 

9 

NC  , 

1. 

1.  to) 

CAwL  LINE 

(CMS  . CyLU 

9 

NC  , 

1. 

1.  4 ) 

CauL  LINE 

(CW2  , CyLL 

9 

NC  , 

1. 

1 . 1 j) 

CAlL  line 

(CW2  , CyLU 

9 

NC  , 

1. 

1.  1 4 ) 

GO  TO  22 

21 

CALL  LINE 

( X . T # l<y . 

i 9 

-1  » 3 

) 

DO  72  L 

= 1 »nd 

B-l 

USL 

No. 
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IF  (GJR(L) .LT.0.5)  GO  10  72 


XL(1) 
XL  (2) 
XL  (3) 
YL(1)= 
YL(2) 


= X ( L ) 

= 0.0 
= 1.0 
Y(l)  U.o 
= tST 
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YL(3)  = 

10.0 

CALL  LINE 

( XL 

* YL 

*i 

*1,1*21) 

72  CONTINUE 

CALL  LINE 

(CSS  , 

Coll 

* NC  , 

1*  1*  26) 

CALL  LINE 

(CSS  , 

colu 

» NC  , 

1#  1#  4 ) 

2«i  CAuL  AXIS 

(0.0  , 

u.o 

#22hPR0PAGATI0N 

>S  IN  Db*22,10.0*9U.O. 

Y (NO  +2)  #10.0  ) 

) 60  TO  15 

(0,0#  u.u#  9HSEA  STATE*9»12.0,0.0,X(N0+1) ,X(NU+2)  *10.0) 


1Y0JD+1)  . 

IFU.Eo.2 
CAlL  axis 
GO  TO  23 

15  CALL  AXIS  (0. 0*0.0. 27HM I NO 
1 XX(N0+1) ,XX(NU+2) ,10.0  ) 

23  CALL  SYMBOL(2.U,9.0,0.l4»lbHPROPA&AUON 
CALL  SYMBOL (3.0*0.5*0.10#  bH VERSUS* 0 . 0 * 6 
IFII.EQ.2  ) GO  TO  14 
CALL  SYMBOL  (2.0,8.0,0.14*  lfltH  SEA 
GO  TO  24 

CALL  SYMBOL (2. 0 ,6.u*0.14* ItiH  WIND 
CALL  SYMBOL (2.0.7.5*0,14* 1&HFREqULNl Y 
IF ( I .EQ.2  ) GO  TO  61 
CALL  SYMB0L(1. 5.7.00 

call  symbol ( 1.5,6.75 

IF  (ISEA.EQ.l)  Go 
CALL  SYMB0L(2. 0,7.0# 
call  Symbol (2. 0,6. 75 

GO  TO  104 

call  symbol (2. 0,7. 00 

1 0.0,  37) 

CALL  SYMBOL (2. 0,6 .75 
1 0.0,  37) 


speed  (mills  per  hour) *27*12. *0., 

LOSS *0.0 #18) 

) 

STATE  *0.0*16) 


14 

24 


SPEED 

1700 


*0.0*16) 
HX*0.0*16  ) 


100 


.10  #4,0. 0,-1) 

.10*26,0.0,-1) 

TO  100 

1.32HPRED1CTI0N  FOR  NEGAT  VE  GRADIENTS. 0*32) 
.1,31HPREOICTIOn  FOR  ISOTHERMAL  WATER #0,0 #31) 

.1,3/HLOWLR  PRLuICTION  FOR  ISOTHERMAL  WATER* 

,1*37HUPPLK  PREDICTION  FOR  ISOTHERMAL  WATER* 


) 


104 

call 

SYMBOL 

) 

CALL 

SYMBOL 

(2,0*6,50#  .10#  13HM£ASUKEl*  LOSS, 0.0,13 

CALL 

SYMBOL 

( 4.0,6.50, ,1*  5,0. 0,-1 

) 

call 

SYMBOL 

( 4.0,6.60, ,1*21*0. 0,-1 

) 

CALL 

symbol 

(4.50,6.50, ,1#12HMINIMuM  LOSS, 0.0*12) 

GO 

TO  69 

ol 

CALL 

SYMBOL 

( 1.5  *7.00*  .10  , 4, 

0.0 

»-l) 

call 

SYMBOL 

( 1.5  *6.75#  .10  , 28, 

u.o 

»-l) 

CAlL 

SYMBOL 

( 1.5  *6.50*  .10  . 14, 

u.o 

»-l> 

CALL 

symbol 

( 1*5  *6,25#  ,10  , 15, 

u.o 

»-l) 

IF 

( isea.eu. 

1)  Go  TO  102 

CALL 

SVMBOH2 

.0*  7.0,.l*47HPREOlcTlON 

T 

c 

X 

NEGATIVE 

GRADIENT (REF tREN 
. 1 , 4ompREUICT lOi.  Foh  ISOTHERMAL  WATLR (REFLRtN 
l.blHPREDlcTIOu  FOR  NlgAUVl  GRADIENT (BLoCK  I 
,1#5UHPREDICT1GN  FOR  lSOfnE«MAu  watER(oLvCK  1 


ICE  ) #0,0*47) 

CmlL  Symbol (2.0, 6. 75 
ICE  ) *0,0* 46  ) 

CAlL  SYMBOL (2. 0,6. 5# 

ISLAND  DATA) #0.0,51) 

CalL  SYMU0L(2.u#6.25 
ISLmNU  DATA) ,0.0,50) 

GO  TO  105 

102  CALL  aYMbOL(2.o, /.00#  .l#52HL0WEr(  PHLuIWION  FO* 
lEFcRENCE  )#  0.0,52) 
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CALL  SYMBOL(2.u,6.75» .1.52huPPtH  HRtulCTlON  FOR  ISO 1 HEKMAl  wAIEH(H 
ltFtRENCE  )>  0.0,52) 

CmlL  SYMdOL(2.o,6.bO« .1»66mlOMLR  PKlUILTIOn  FOK  ISO rhLuMAL  wArcR(o 
1LOCK  ISLAND  DATA) »u. 0,56) 

Call  SYMBOL (2. 0,6 *26*  «l»56FtuPPEh  PWLulcriON  FoK  ISOTHERMAL  WA  ILH  (b 
1LOCK  ISLAND  DATA ) » u . 0 » 56) 

106  CALL  SYMBOL  ( 1.5,7.25, ,10, 5,0. 0,-1  ) 

Call  SYMBOL  (2.0» /.25»  ,10,13HMEASUHcU  LOSS. o.o, 13  ) 

CAuL  SYMBOL  ( 4.0,7.26. .1.  5,0. 0,-1  ) 

CALL  SYMBOL  ( 4,0,7.35, , 1,21,0. 0,-1  ) 

CALL  SYMBOL  (4 , 50 , 7.25, , 1 , 12HM1N1MUM  LOSS, 0.0*12) 
o9  IF ( I ,EQ,2  ) GO  TO  20 
DO  7 K=  1,  NS 
FN(K)  = 0.0 
FM(K)  = 0.0 

V CK ) = 0.0 

w* IK)  =0.0 

MMlN  ( K ) = 0.0 

WWwl(K)  = 0.0 
WWw2(K)  = 0,0 


GTA(K) 

= 0 

. o 

7 

W (K)  = 

0.0 

GO  TO 

26 

20 

DO  47 

K=  1, 

NW 

FN(K)  = 

0.0 

FM  (K ) 

= 0.0 

V (K)  = 

0,0 

WW  (K) 

=0.0 

WWW (K) 

= 0. 

0 

WWWl (K) 

= 0 

.0 

WWW2(K) 

= 

0.0 

GTA(K) 

= 0 

.c 

47 

W (K)  = 

0.0 

DO  4 J 

si,  NO 

IF  ( I ,£<i 

.2  ) 

GO  TO  16 

K = 2.0 

*X(J) 

+ 1.0 

WHITE 

( 4,95 

> K 

GO  TO 

25 

16 

XX(J) 

= XX(J>/  5,0 

WHITE  (4,96)  AX(7) 

WHITE  (4,96)  XX(6, 

K =2.Q*XX(J)  +1,0 

WHITE  ( 4,95  ) K 
26  IF (K.GT .3)  K=  K - k/3 
WHITE  ( h,95  ) K 
IF (K.GT .6)  K=  A — i 
WHITE  ( 4,95  ) K 
IF(K,GT.10)K=  K - i 
WHITE  ( 4,95  ) a 
95  FOHMAT  (IS) 

F.M(K)  = Fn(K)  + 1.0 

IF  (GTN(j) .OT.0,5)  GTA(a)  = I.o 
w(K)  = w(K)  + Y(J) 

Ww  (k ) = WW(K)  4 YY ( J) 

V(K)  = V(K)  ♦ Y (0)4*2. 0 
IF  (GTR( J) .G1 ,0.3)  60  TO  99 

FM(K)  = FM(K)  > l.o  USL  Tech  Memo 

w*a(K)  = WWW (K ) 4 YYT(J)  r , No.  2211-91-69 
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*1*11*1  (K ) = wWWl(k)  ♦ YYYKK) 

' MMV.2  (K ) = *IW«.2(K)  ♦ YYY2(K) 

'J'J  *r<lTE  (4,11)  FN(K)  * W(K)  , Y(K)  » K 

11  FOKMAT  ( 3F10.1  , 110  ) 

4 CONTINUE 

IF ( I ,EU,2  ) GO  TO  ltt 
READ  (3,6) (2(K)  , K=  1*NS  ) 

GO  TO  2o 

lo  READ  (3,6)(ZZ(K),  K=  1 * N'W  ) 

Zo  IF ( I ,LQ,2  ) GO  TO  17 
DO  5 K=  1*  NS 

IF  (FN(K) .OT.0.5)  GO  TO  ZOO 
FN(K)  = FN(K-l) 

WOO  = W(K-l) 

WW(K)  = W*(K-1) 

Z ( k ) = Z(K-l) 

V(K)  - V(K“1) 

GO  TO  5 

200  WOO  = W(K)/FN(K) 

W*(K)  = WW(K)/FN(K) 

5 S(K)  = SORT ( V(K)/FN(K)  - W(K)**2.0  ) 

GO  TO  60 

17  DO  45  K=  1*  Nw 

IF  (FN(K) .GT.O.b)  GO  TO  201 
FN(K)  = FN(K-l) 

W(K)  = W (K“l ) 

WW (K)  = WW(K-l) 

ZZ(K)  = ZZ(K-l) 

V (K)  = V(K-l) 

GO  TO  45 

201  W (K)  = WOO/FNOO 
WW(K)  = WW(K)/FN(K) 

45  S(K)  = SORT ( V(K)/FN(K)  - »(K)**2.0  ) 

6o  call  plot  (20.0,0.0,  -3) 

6 format (8fio,5) 


IF ( I ,EO, 

2 ) GO 

TO  19 

DO  12 

K = 1* 

NS 

SI (K ) = 

W (K) 

+ S(K) 

12 

S2(K)  = 

W(K) 

- S(K> 

GO  TO 

27 

IV 

DO  42 

H 

11 

* 

Nft 

S1(K)  = 

w(K) 

♦ S(K) 

42 

S2(K)  = 

W (K) 

- S(K) 

IF (l.EO, 

2 ) GO 

TO  62 

27 

W ( NS+1 ) 

= ESI 

W (WS+2) 

=10.0 

SKNS+1) 

= EST 

S11NS+2) 

= 10.0 

bZ (NS+l ) 

= EST 

S2(NS+2) 

= 10.0 

CS5<nC+1)  = 0.0 
CSS ( NC+2 ) = 1.0 
Z(inS+1)  = O.U 
Z(NS-fZ)  = 1.0 


CALL  LINE  (Z*WMS*  1,1,  5 ) 

DO  81  l=  1 * NS 

if  (GTA(L).LT.0.5)  oo  ro  bi 

USL 

Tech  Memo 

ZL(1)=  Z(L) 

No. 

2211-91-69 
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ZL(2)  = U.O 

ZL(3)  = 1.0 

WL(1)=  *IL>  ♦1.0 
*L(2)  = EST 

*1(3)  = 1U,0 

call  line  ( ZL  , WL  ,1  ,1,1,21) 

61  CONTINUE 

CALL  LINE  (Z.Sl.NS,  1,-1,  22) 

CALL  LINE  (Z.S2.NS,  1,-1,  23) 

CALL  LINE  (CSS»COLu,NC,  1,  1»  4 ) 

CALL  LINt  (CSS.COLl.NC,  1»  1»  26) 

CALL  AXIS  (0.0»  0.0.  9HSEA  sTATE»9 » 12. 0 , 0 . U . X (ND+1 ) ,X (N042)  .10.0) 
CALL  AXIS  (0.0  , u.O  »22rtPR0PAGATION  LOSS  IN  GB» 22 , 10 . 0 » 90 . 0 , 
lY(NO-fl)  . r ( ND  42)  .10.0  ) 

CALL  SYMBOL(2.O,9.u,O.14»10hPR(JPA6ATION  LOSS. 0.0. 16) 

CALL  SYMBOL(3.0,B.5,0.10»  6HVEkSUS»0.U,b  ) 

CALL  SYMBOL (2. 0.0, u.O. 14. 16H  SEA  STATE  .0.0.18) 

CALL  SYMBOL (2. 0.7. S. 0.14. 16hFREuUEN(,  Y 1700  HZ»Q.0,16  ) 

Call  SYMBOL(1.S,7.uO».10  ,4,0. 0,-1) 

CALL  SYMBOL (1.5,6.75, .10,26,0. 0,-1) 

IF  (ISEA.EQ.1)  Go  TO  101 

CALL  SYMBOL(2.O,7.0.  ,l»32hPREDIcTIO|.  FOR  NEGATIVE  oRAD1anT#0.G»32) 
CALL  SYMBOL (2. G.b.  75.  ,1. 3lHpREDIC  T ION  F JR  ISOTHERMAL  WATER. 0 • 0 . 31 ) 
GO  TO  106 

101  CALL  SYMBOL (2. 0,7.00. «1,37HlOWER  PREDICTION  FOR  ISOTHERMAL  WATER. 

1 0.0,  37) 

CALL  SYMBOL(2.0,b.7S» .1.37MUPPER  PREulCTlON  FOR  ISOTHERMAL  wATtR, 

1 0.0,  37) 

10b  CAlL  SYMBOL  (1.5.6,50. .10.5,0. 0,-1) 

CALL  SYMbOL  (2.0.6,50. ,10. l4HAVERAGt  VALUES,  0,0.14) 

CALL  SYMbOL  ( 4.0,6.50, ,1.  5,0. 0,-1  ) 

CALL  SYMbOL  ( 4. 0 , 6,60 , , 1 » 2l » 0 . 0 , -1  ) 

CALL  SYMBOL  ( 4 , 50 , 6.50 , , 1 , 12HMINIMUK  LQ5S»0.0,12) 

CALL  SYMbOL  (1.5,6.25, .10,22,0. 0.-1) 

CALL  SYMBOL  (2 . 0 » 6 . 25.  . 10  » 2<4HUPPEk  STANDARD  DEVIATIONS. 0,24) 

CALL  SYMBOL  (1.5,6.00, . 10»23»0.0»-1) 

CaLL  symbol  (2.0,6.00, .1Q»24HL0WER  standard  DEVIATION, 0.0, 24) 

WRITE  (4,10) 

10  FORMAT  ( 50H  StA  SiATE  AVER  LOSS  NUM  VAlVES  STANDARD  DEVIATION) 

DO  6 M=  l , NS 
9 format  <4Fio.5) 

a WRITE  (4,9)  Z(N')  » W (M)  , FN(M)  , S(f) 

IF  (ISZ.EO.l)  Go  TO  30 
W* ( NS  +1)  = EST 

WW ( NS  +2)  s 10. 0 
( NS  41)  = 0.0 

WWw(NS  42)  = 10. u 

WWtaltNS  41)  = G.u 

*W»1(NS  42)  = 5.0 

www2 ( NS  41)  = 0 • o 

*w„2(NS  42)  = 5,0 


CALL  PLOT  (20.0,  u.O  » -3  ) 
CALL  LINE  (Z  » W,  NS, 
CALL  LINE  (Z  , WW,  NS, 

DO  o2  L=  1 » NS 
IF  (GTA(l) .LT.0.5)  GO 
ZL(1)=  / ( L ) 

ZL ( 2 ) = U.O 


1 

1 


,1 

» 1 


3 ) 

4 ) 


TO  62 
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ZL  ( 3 ) = 1.0 

mL(15=  M Cl)  +1.0 
ML (2 ) = LST 

ML (3)  = 10,0 

CALL  LINE  ( ZL  * ML  *1  *1,1. 2l) 

02  CONTINUE 

CALL  AXIS  (0.0  , u.O  * 22HPR0PAGAT ION  LOSS  IN  08*22,10.0*90.0, 
1Y(ND+1)  , Y(ND  +2)  ,10.0  ) 

CALL  AXIS  (0.0,  O.o*  9HSEA  STATt * 9 , 12 . 0 , u . 0 , X ( ND+1 ) , X ( l<L)+2 ) *10.0) 
CALL  SYMByL (2.0,9.0,0.14* 18HPR0PAGAT ION  L05S»0.0*16) 

call  Symbol (3. o,8. b*o, io»  6hveksus* 0.0*0  ) 

CALL  SYMBOL(2.0,8.0.0.l4*lbH  SEA  STATL  *0.0*18) 

CALL  SYMBOL(2.0,7.S.0.14»18HFKLQUENCY  1700  HZ»0.0*18  ) 
call  Symbol (2. 0,7 .oo* .10  *4,o.o,-i> 

CALL  SYMBOL ( 2. s, 7.0 * ,1*32HMaXIM0M  OvEK  PULSE  *0.0*32) 

call  symbol (2. o,b.75*. 10*  b,o.o,-i) 

call  symbol (2. s, 6. 75, ,i*3ihaveka&l  over  pulse  *o.o*3d 

CALL  symbol  (4.75,6.75, ,1,  5,0.0,-l  ) 

CALL  SYMBOL  (4.75,6.85, ,1*21*0. 0,-1  ) 

CALL  SYMBOL  (5.00,6.75,,1,12HMINIMUM  LOSS»0.0»12) 

CALL  PLOT  (20.0*  0.0  * -3  ) 

DO  205  K=  1 * NS 

IF  (FM(K) .GT.O.b)  60  TO  202 

NMM(K)  = MMW(K-l) 

MWMl(K)  s mWWI(K-I) 

WMM2 ( K ) = MWW2 ( K-l ) 

Z(K)  = Z(K-l) 

60  TO  205 

202  WMM(K)  s MMW(K)/FM(K) 

MMmI(K)  = WWM1(K)/FM(K) 

MWM2 ( K ) = MWM2(K)/FM(K) 

20b  CONTINUE 

CALL  LlNt  (Z  , mWM * Na,  1 *1  , 4 ) 

CALL  AXIS  (0. 0*0.0,  8HVaRIAnCE,6* 10 . ,90 . , MMM (NS+l ) , *MM (NS+2) » 10,0) 
CALL  AXIS  (0.0*  o.o*  9USEA  STATE  * 9* 12. 0 , 0 . 0 , X (ND+1 ). X (NO+2)  *10.0) 
CALL  SYM80L(2.0,9.0,0.14*18h  VARIANCE  *0.0*18) 

CALL  SYMBOL (3. 0,6. 5*0.10*  6HVEKSUS* 0.0*6  ) 

CALL  SYMBOL(2.0,8.0*0.14*18h  SEA  STATt  *0.0*iB) 

CALL  SYMBOL(2.0,7.5*0.14*18hFRLQUENi.Y  1700  HZ*0.0»18  ) 

call  symbol (2.0,7.00**10  *4,o.o,-i) 

CALL  SYMB0L(2.b,7.u* ,1*32HVaRIANCE  OVtK  PULSE  *0.0*32) 

CALL  PLOT  (20.0,  o.O  * -3  ) 

CALL  LINE  (Z  , WmMI*  NS*  1 *1  , 5 ) 

CALL  LINE  (Z  , M*M2*  NS*  1 *1  , 4 ) 

CALL  AXIS  (0. 0*0.0,  8HVARIAnCE,8*10. ,90. ,mmm(NS+1) *„MM(NS+2> *1U,o) 
CALL  AXIS  (0.0*  0.0,  9HSEA  STATE *9* 12.U , 0 .0 , X (NU+l > »X (ND+2)  *10,0) 
CALL  SYMBOL (2 . 0, 9.0*0.14, 18h  VaKI ANCt  *0.0*16) 

call  symbol (3. 0,8.5*0.10*  6hversus*u.o,o  ) 

CALL  SYMBOL (2. 0 ,8.0*0. 14*1 6h  SEA  STATt  *0.0*18) 

CALL  SYMBOL(2.0,7.5,0,l4*lttHFREQULNLY  1700  HZ*0.0*18  ) 

CALL  SYMBOL (2. 0 ,7.00*.l0  *4,0. 0,-1) 

CALL  SYMHOL(2.b,7.o* ,l*32HVARIANCt  OF  AVtKAoE  OVER  PULSE  *0.0*32) 

call  symbol (2. 0,6. 75* . 10*  s,o.o,-d 

CALL  SYMBOL (2. 5,6, 75, . 1 , 31HVARI ANCE  OF  MAXIMUM  OVLK  PULSL  ,0.u*3l) 
60  TO  30 

62  ZZ(Um  ♦ 1)  5 0.0 

ZZ(NW  + 2)  = b.o  USL  Tech  Memo 

W(NM  ♦!)  = EST  B-6  No.  P211-91-69 
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M(hM  ■ 

4-2)  = 10.0 

S1(NM 

♦1)  = EST 

SI  (NW 

♦2 ) - 10.0 

S2(NW 

+1)  = EST 

S21NW 

■•■2)  = 10.0 

CWStNC 

+ 1)  = 0.0 

CMS (NC 

♦2)  = 5.0 

CM2 (NC 

♦1)  = 0.0 

CM2 (NC 

4-2)  = 5.0 

call  1 

LINE  (ZZ  , M ,NM  , 

1. 

1.5) 

DO  83 

L=  1 » NW 

IF  (GTA(L) .LT.0,5)  GO 

TO 

83 

ZZL(l) 

= ZZ (L) 

ZZL ' 2 ) 

= o.o 

ZZL(3) 

= 5.0 

ML(1)= 

W(L)  +1.0 

ML  (2) 

= EST 

ML  ( 3) 

= 10.0 

CALL  LINE  (ZZL  . ML 

.1 

*1 

USL  Tach  Memo 
No.  2211-91-69 


83  CONTINUE 
CALL  LINE 
CALL 
CALL 
CALL 
CALL 
CALL 


»l.l*2l) 


LINE 

LINE 

LINE 

LINE 

LINE 


»Si 

*S2 


NW  .1,-1. 22) 
NW  ,1,-1, 23) 


(ZZ 
(ZZ 

(CMS  » COLU  » NC  .1,1.4) 

(CMS  . COLL  . NC  .1,1.28) 

(CW2  r COLU  » NC  .1,1,14) 

(CM2  ' COLL  t NC  ,1,1,1b) 

CALL  AXIS  (0.0'0.0,27HWINO  SPEED  (MILES  PEK  HOUR) #27, 12. ,0. , 

1 XX(ND+1) ,XX(ND+2),10.0  ) 

CALL  AXIS  (0.0  , 0.0  .22HPR0PAGATI0N  LOSS  IN  DB.22, 10.0,90.0, 

1Y(ND+1)  , Y(NO  +2)  .10.0  ) 

CALL  SYMOoL(2.U.9.0.0.14»18hPROPAGAT1ON  LOSS, 0.0, 18) 

CALL  SYMBOL (3. 0,8.5,0.10)  6HVERSUS,u.0,b  ) 

CALL  SYMBOL(2.0,b.0,0.14»18H  WIND  SPEED  ,0.0,18) 

CALL  SYMB0L(2.0,7.5,0.14,18hFRE&UENlY  1700  HZ’0.0»18  ) 


CALL 

SYMBOL 

( 

l.b 

.7.00. 

.10  , 

4, 

0.0 

.-1) 

CALL 

SYMBOL 

( 

1*5 

,6.75. 

.10  . 

28, 

0.0 

»-l) 

CALL 

SYMBOL 

( 

1.5 

.6.50. 

.10  . 

14, 

0.0 

*-l) 

CALL 

SYMBOL 

( 

l.b 

.6*25. 

.10  , 

15, 

0.0 

.-1) 

IF 

(ISEA.EQ. 

1) 

60 

TO  103 

CALL  SYMBOL(2.0,7.0,.1»47HPRED1CTIOn  FOR  NEGATIVE  GRADIENT (REFEREN 
ICE  ) >0.0.47) 

CALL  SYMBOL (2.0.6. 75 » , 1»4ohpREUICTIun  FOR  ISOTHERMAL  MATEH(REFEREN 
ICE  ), 0.0,46  ) 

CALL  SYMBOL(2.0,b.b».l*5lHPREDlCTXON  FOR  NEGATIVE  GRADIENT (BLOCK  I 
ISLAND  DATA) ,0.0,51) 

CAlL  SYMBOL (2. 0,6. 25. .1,50HPREUICTIDn  FOR  ISOTHERMAL  MATEK(BLoCK  I 
ISLAND  DATA) .0.0, SO) 

GO  TO  107 

103  CALL  SYMBOL(2.0,7.U0..1.52HLOWER  PREDICTION  FuH  ISOTHERMAL  WA TER (R 
1EFERENCE  ).  0.0,52) 

CALL  SYMBOL (2. 0,6. 75. . 1.5ZHuPPER  PRldICTION  FuR  ISOIHEkMAl  MATER(R 
1EFERENCE  ).  0.0,52) 

CALL  SYMBOL (2.0, 6. b0..1,5oHL0WER  PREDICTION  FOR  ISO  THERMAL  MATER (B 
1LOCK  ISLAND  DATA) .0.0.56) 

CALL  SYMB0L(2. 0,6.25.  ,1»56HuPPER  PREDICTION  for  ISOJHlkMAl  wAlu.R(U 
1LOCK  ISLAND  DATA) »u, 0.56) 

107  CALL  SYMBOL  (1.5. 6, 00. .10.23*0.0.-1) 

CALL  SYMBOL  (2.0.b,00. .10.24HLOMER  STANDARD  DlV  I A UuN. U .0 ,24 ) 

B-7 
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CALL  SYMBOL  ( 1 . 5# 5. 75# . 10 # 22* 0. 0 *-l ) 

CALL  SYMBOL  (2.0*5,75* ,10*24HUPPER  STANDARD  DEVIATIONS. 0,24) 
CALL  SYMBOL  (1.5,5.50, .10, 5,0. 0,-1) 

CALL  SYMBOL  (2.0»5.50# . 10 # l4HAVERAGt  VALUES.  0.0#l4) 

CALL  SYMBOL  (4,25.5.55. ,1#  5.0. 0,-1  ) 

CALL  SYMBOL  (4.25,5.65, (l»2l,0. 0,-1  ) 

CALL  SYMBOL  (4.50,5.55, (1,i2HMIN1M0M  LOSS#0.0*12) 

WRITE  (4,90) 

90  FORMAT  ( 50H  WIND  SPFED'VER  LOSS  NUM  VALVES  STANDARD  DEVIATION) 
DO  80M=  1 , NW 

BO  WRITE  (4,9)ZZ(M)  , W(M),  FN(M)  , S(M) 

IF  ( ISZ.EU, 1 ) Go  TO  JO 

WW (NW  +1)  = EST 

WW (NW  42)  = 10. 0 

WWW ( NW  +1)  = 0.0 

WWW (NW  +2)  = 10.0 

WwWl (NW  +1)  = 0.0 

WWW1(NW  42)  = 5.0 


WWW2(NW  -t-1)  = 0.0 

WWW2(NW  42)  = 5.0 

CALL  PLOT  (20.0#  0.0  # -3  ) 

CALL  LINE  (ZZ,  W,  NW#  1 #1  , 5 ) 

CALL  LINE  (ZZ,  WW#  NW#  1 #1  # 4 ) 

DO  84  L=  1#NW 

IF  (GTA(L) .LT.0,5)  GO  TO  84 

ZZL(1)=  ZZ(L) 

ZZL (2)  s 0.0 
ZZL(3)  = 5.0 

WL(1)=  W (l)  +1.0 

WL (2)  = EST 

WL13)  = 10,0 

CALL  LINE  (ZZL  # WL  #1  #l,l#2l) 

84  CONTINUE 

CALL  AXIS  (O.o  , O.O  #22HPr0PAGATI0N  LOSS  IN  DB» 22, 10,0# 90.0, 
lY(ND4l)  , Y (ND  42)  #10.0  ) 

CALL  AXIS  (0, 0*0.0, 27HW XNO  SPEED  (MILES  PER  HOUR) ,27*12. ,0. . 

1 XX(ND4l) ,XX(ND+2) ,10.0  ) 

CALL  SYMBOL(2.0,9.0#0.14»18hPROPAGATION  LOSS#0.0*l8) 

CALL  SYMBOL (3.0,8.5,0.10,  OHVERSUS# 0.0*6  ) 

CALL  SYMBOL(2.0, 8.0»0. 14, 18H  WIND  SPEED  *0.0*18) 

CALL  SYMBOL (2.0, 7.5*0.14, IBhFREQUENc Y 1700  HZ*0.0*18  ) 

call  symbol (2. o,7. oo* .10  *4,o.o,-n 

CALL  SYMB0L(2.5,7.U,.1*32HMaXIMUM  OVER  PULSE  *0.0*32) 

call  symbol (2. o,6 .75* .10,  5,o.o,-d 

CALL  SYMB0L(2.5,6.75*.1,31HaVERaGE  OVER  PULSE  *0.0*31) 

CALL  SYMBOL  ( 4.5,6.75, ,1*  5,0. 0,-1  ) 

CALL  SYMBOL  ( 4,5,6.85, ,1,21,0. 0,-1  ) 

CALL  SYMBOL  (5, 00 ,6.75, , 1 , 12HMINIMUM  LOSS, 0.0*12) 

CALL  PLOT  (20.0*  0.0  * -3  ) 

DO  245  K=  1*N* 

IF  (FM(K) .GT.0.5)  GO  To  203 
WWW (K ) = WWW(K-I) 

WWwl(K)  = WWWl(K-l) 

WWW2 (K ) = WWW2 (K-l ) 

ZZ (K ) = ZZ (K-l ) 

GO  TO  245 

203  WWW(K)  = WWW(K)/FM(K)  USL  Tech  Memo 
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MMM2(K)  = *W*2(K)/FM(K) 

24b  CONTINUE 

CALL  LINE  (ZZ,  ftMM,  NM,  1 *1  , 4) 

CALL  AXIS  (Q.O.O,O,0HVARXANcE,8,lO,,9O,,MMMl(N«i4l)»*MMl<NM+2> ,10.) 
CALL  AXIS  (0.0»0.0,27HMIND  SPEED  (MIlES  PER  HOUR) ,27.12. ,0. * 

1 XX(ND+1) ,XX(ND+2) ,10.0  ) 

CALL  SYMBOL(2.0,9.0.0.14»1BH  VARIANCE  »0.0»i8) 

CALL  SYMBOL(3.0,8.5»0.l0»  6HVERSUS. 0 . 0 , 6 ) 

CALL  SYMBOL(2.0,8.0'0.14»16H  MIND  SPEED  '0.0*16) 

CALL  SYMbOL(2.0.7.b'0.14»18HFREQUENCY  1700  HZ'0.0'18  ) 

CALL  SYMBOL (2. 0.7* 00' .10  .4,0, 0,-1) 

CALL  SYMBOL (2. 5' 7*0' . 1.32HVARIANCE  OVER  PULSE  '0.0 '32) 

CALL  PLOT  (20.0*  0.0  , -3  ) 

CALL  LINE  (ZZ.  MMM1,  NM,  1 »1  , 5 ) 

CALL  LINE  (ZZ.  MftM2,  NM.  1 .1  , 4 ) 

CALL  AXIS  (0.0'0.0.8HVARXANcE'd.l0*.90..MMMl(NM+l) ,MMMl(NM+2) .10.) 
CALL  AXIS  (O.O'O.O. 27HM I ND  SPEED  (MILES  PER  HOUR) .27' 12. . 0. ' 

1 XX(NO+l) ,XX(ND+2) ,10.0  ) 

CALL  SYMBOL (2.0. 9.0. 0.14. 18h  VARIANCE  .0.0.18) 

CALL  SYMBOL (3*0, 8.b. 0*10.  6hVERsUS»0,0.6  ) 

CALL  SYMBOL(2.0,8.0.0.14.18h  MIND  SPEED  .0.0.16) 

CALL  SYMOOL(2.0,7.5»0.14»18HFREOUENCY  1700  HZ' 0.0 .18  ) 

CALL  SYMBOL (2. 0,7 .00. .10  »4,0.0,-l) 

CALL  SYMBOL(2«5,7.0,.1'32HVaPIANCE  OF  AVERAGE  OVER  PULSE  .0*0.32) 
CALL  SYMBOL (2*0,6*75. .10,  5.0.0.-1) 

CALL  SYMB0L(2*b,6.75..1'31HvARlANCE  OF  MAXIMUM  OVER  PULSE  .0.0*31) 
30  CONTINUE 

CALL  PLOT  (20.0.0.0.  -3) 

CALL  PLOT  (0.0.0.0,900) 

ENO 
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